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ABSTRACT 



Aims. Intermediate Mass (IM) stars are an important component of our Galaxy, as they significantly contribute to the interstellar 
FUV field and, consequently, play an important role in the energy balance of the ISM. Despite their importance, very little is known 
about their formation process and only a few studies have been devoted to characterize the first phases in the evolution of intermediate 
mass protostars. Here we consider in great detail the case of the brightest and closest known young IM protostar: FIR4 in the 0MC2 
component of the Orion molecular cloud complex. 

Methods. We analyzed the available continuum emission (maps and SED) through one-dimensional dust radiative transfer calcula- 
tions. We ran large grids of models to find the envelope model that best fits the data. The derived dust density and temperature profiles 
have been then used to compute the gas temperature profile, equating gas cooling and heating terms across the envelope. Last, we 
computed the water line spectrum for various possible values of water abundance. 

Results. The luminosity of FIR4 has been reevaluated to 1000 L©, making F1R4 definitively an Intermediate Mass protostar The 
envelope surrounding FIR4 has a relatively shallow density power law index, ~ 0.6. The most surprising result is that the gas and 
dust are thermally decoupled in the interior of the envelope, where the dust ices sublimate at 100 K. This has important consequences 
in the interpretation of the line data. We provide the predictions for the water spectrum, and discuss in detail the lines which will be 
observed by the Herschel Space Observatory . 

Key words. ISM: abundances — ISM: molecules — stars: formation 



1. Introduction 

, Intermediate mass (IM) stars, namely stars whose mass is in 
the 2 to 8 Mq range, are crucial in studies of star forma- 
tion becau se they constitute the Unk between low- and high- 
mass stars jPi Francesco et alJl997tlMannings & Sargentll997l 
12000) . and, therefore, can help to understand if and how much 

I different are the processes at work in the two ends. On the one 
hand, low mass stars are can be formed isolated or in loose 
groups of few objects per cubic parsec (Gomez etal. 1993), 
while hi gh-mass stars are usually found to form in tight clus- 

' ters (e.g iHillenbrand & Hartmannlll998l). IM stars, on the other 

' hand, are also found in cltisters (e.g. iTesti & SargentI Il998t 
iNeri et alj2007l:lFuente et alj2007h . with a smooth transition to- 
wards the low mass star, loose c l uster regim e for star masses 
around 3.5 Mq dTesti et al.lll999l) . iTesti et all (l99% also con- 
cluded that IM stars mark the transition from low density aggre- 
gates of < 10 stars per cubic parsec of T Tauri stars to dense 
clusters of > lO"* stars per cubic parsec associated with early- 
type stars. In agreement with the different observed environ- 
ments, several authors have proposed that high mass stars are 
formed by coalescence of lower mass stars, whereas other au- 
thors favor the "monolithic" f o rmatio n (see for example the re- 
cent review by iBeuther et alj (l2007h ). In this context, the IM 
stars study can greatly help the debate. Indeed, due to their 
intermediate position, the study of IM protostars will provide 
crucial information on the transition between the two forma- 
tion regimes as well as on the limits of the low mass and high 
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mass formation scenarios. Finally, IM stars are ar nong the dom- 
inant sources of the Int er-Ste llar FUV field (e.g. iHabindl 1 9681 : 
Gondhalekar & Wilson 19751) , which regulates the phases of the 
ISM in the Galaxy, and, in turns, the overall Galaxy star forma- 
tion process and history. Despite the far-reaching importance of 
IM stars, very little is known about the formation and first evolu- 
tionary stages of these stars. The situation is so bad that to date 
we do not have a satisfying sample of Class IM objects, namely 
objects representing the first phases of stellar formation, where 
the protostar is embedded in its envelope and its luminosity is 
dominated by the accretion luminosity, nor a systematic study 
of their physical structure, as it is the case for low mass Class 
sources (e.g. Ceccarelli et al. 2007; Di France sco et al . 2007j). 
This article is the first of a series that aims to fill this gap in our 
knowledge. 

In this context, the Orion Molecular Cloud 2 (OMC2), the 
closest known region where high to low mass star formation 
is going on, represent s a precious laboratory for these stud- 
ies. Observed first by iGatlev et all ( fl974.) , OMC2 is located 
15' (~ 2 pc) North of the Orion nebula. It has a filamentary 
structure, elongated in the direction north-south, with active 
star formation concentrated in the central and densest region, 
shielded from the U V radiation from newly formed OB stars 
dJohnsoiLet al.Hl990^he mass of the cloud amounts to about 
1500 Mq (Mezger etal]|l990l) . Several extensive studies have 
shown that OMC2 is a rich star forming region, which har- 
bors several young protostars, including sever al Class can- 
didates (A li & Depov 1995; Chinietal. 1997; iLiset all 119981; 
iJohnstone & Bally, 1 999.; .Reipurth et al.i,1999) . Observations of 
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molecular lines have revealed several outflows emanating from 
the young protostars in the region. Many studies have focused 
on the out flows (e.g. [Wi lliams et al. 2003) and their impact on 
the cloud (lAso et al.ll2000t IWu et al.1 12005 ). Only few of these 
studies, in contrast, have addressed the problem of the chemic al 
structure of the forming stars in 0MC2 (iJohnstone et al.ll2003h . 

Among the several protostars i n OMC2, FIR 4 stands out as 
the brightest submillimeter source dMezger et al.. 1990 ). Located 
almost at the center of the clo ud, FIR4 is also a br ight IRAS 
source and a VLA radio source d Reipurth et alJ[T999 l). All these 
characteristics led i Reipurth et al . ( 1999) to define FIR4 "a bona 
fide Class source". The FIR4 integrated luminosity was esti- 
mated to be about 400 Lq and the envelope mass is about 35 
Mq. Such values led to identify FIR4 as an intermediate mass 
protostar (I Johnstone et al.ll2003l) . Because of its vicinity and its 
relatively bright molecular lines, FIR4 is an ideal source for a de- 
tailed study of the physical and chemical structure of an IM pro- 
tostar. Existing dust continuum and molecular line observations 
point to an envelope with at least two components: a warm com- 
ponent with an average temperature of about 40 K and a colder 
component at about 1 5 K (Mezgeretal. 1990; Johnston e et alj 
l2OO3b . lJ0rgensen et alj (2006) modeled the 850 yum SCUBA map 
towards this source to reconstruct its temperature and density 
profiles. Based on the ob s erved CO and H2CO millimeter line 
emission, J0rgensen et al. (2006) concluded that the FIR4 enve- 
lope is illuminated by an external FUV field amounting to 1 x 10"* 
times the Interstellar FUV field. However, their interpretation 
suffers of some contradictions emphasized by the same authors. 
For example, such an intense FUV field would heat up the whole 
envelope to a temperature larger than 25 K, the CO freezing 
temperature (Qberg et aL 2005), in contradiction with the mea- 
sured average CO abundance, ten times lower than the canonical 
value, which would rather testify for a large CO-frozen region 
(lJ0rgensen et a l. 2006). In addition, the maps of the fine struc- 
tur e lines of the O and C^ atoms together with the CO 1-0 line 
led [Herrmann et af] (1 19971) to conclude that the OMC2 region is 
illuminated by a FUV field 500 times the Interstellar field. 

Given this puzzling situation, we decided to derive again the 
temperature and density profiles of FIR4 by taking into account 
more data than those considered bv lJ0rgensen et al.l (2006) (©. 
Using the derived dust temperature and density profiles, we then 
computed the gas temperature profile, by equating the heating 
and cooling term s across the env elope ( @. As sho wn by sev- 
eral authors (e.g.lC eccarelli et al., 1996; Doty & N eufeld 1997), 
the gas cooling in protostellar envelopes is dominated by the 
emission from the rotational lines of CO and, more important, 
H2O together with the fine structure lines of OI. Actually, wa- 
ter is a key molecule in the gas thermal balance for two rea- 
sons. First, in the warm regions where the grain mantles subli- 
mate, it is the most abundant molecule; second, given its rela- 
tively large dipole moment, water is a very powerful line emit- 
ter and, consequently, gas coolant. Given its major role in the 
prediction of the gas temperature profile, we discuss the depen- 
dence of the derived gas temperature on the assumed water abun- 
dance profile, which is poorly known. Not surprising, FIR4 is in 
fact one of the few sources where the full spectrum between 50 
and 2000 GHz is planned to be observed at high spectral resolu- 
tion with the Heterodyne Instrument for the Far Infrared (HIFI) 
on board Herschel (http://herschel.esac.esa.int/), to be launched 
in 2009. HSO, and specifically the high resolution interferom- 
eter HIFI, will allow to observe the water lines in the 500 to 
2000 GHz range with unprecedented spectral and spatial resolu- 
tion. Motivated by the Herschel mission, we report the predicted 
water line spectrum for the different assumed water abundance 



profiles, and discuss the observability by HIFI and PACS (® . 
Section|5]concludes the article. 



2. Dust density and temperature profiles 

In this section, we derive the dust density and temperature pro- 
files by modelling the 350, 450 and 850 fj.m maps of the region, 
plus the Spectral Energy Distribution (SED) from the millime- 
ter to the Mid-Infrared (MIR) wavelength range. We first de- 
scribe the observations we used in our analysis 02.1b and then 
the modeling ( 92.21 ) and the result of the modeling ( 92.3l l. 

2.1. Continuum emission: observational data 

In our analysis, we use the maps of the continuum emission at 
850, 450 and 350 yum obtained at JCMT and CSO respectively. In 
addition, we take into account the Spectral Energy Distribution 
(SED) of FIR4 from 24 to 850 i-im obtained considering also the 
IRAS and Spitzer observations. 

a) 850, 450 and 350 jim maps 

We retrieved the 450 and 850 yum maps obtained by 
IJohnstone & BaUvl 11999.) at the 15 m James Clerk Maxwell 
Telescope (JCMT) with the focal-plane instrument SCUBA 
(Submillimeter Common-User Bolometer Array). The spatial 
resolution of the maps is 7.5" and 14.8" at 450 and 850 pm re- 
spectively. The calibration uncertainty and noise levels are esti- 
mated by those authors < 10% and 0.04 Jy beam"' at 850 pm and 
< 30% and 0.3 Jy bea m"' at 450 yum, respectively. The 350 yum 
map was obtained by Lis et al.l (Il998l) at the 10.4 m telescope 
of the Caltech Submillimeter Observatory (CSO). The instru- 
ment used was the bolometer camera SHARC. The resolution of 
the map is 12". The calibration uncertainty has been evaluated 
~ 25%-30%. The three maps are reported in Fig.[T] They show 
the envelope surrounding/forming FIR4 which extends for about 
20", but also the presence of two sources: FIR3, 25" North, and 
FIR5, 25" South. To evaluate the continuum brightness profile 
of the FIR4 envelope, we averaged the continuum flux over ring 
at the same distance from the FIR4 center, excluding the regions 
contaminated by the presence of FIR3 and FIR5 (dashed regions 
in Fig.[T]). The resulting brightness profiles are shown in Fig.|2] 
Note that in the analysis of the envelope emission 02.2t we sub- 
tracted the cloud contribution, estimated to be ~ 0.001, ~ 0.03 
and ~ 0.05 Jy arcsec"^ at 850, 450 and 350 yum respectively. 
Furthermore, in order to take into account that the SCUBA and 
SHARC maps were obtained with the chop throw of 65" and 
90-120" respectively, we only considered the inner 60" in our 
analysis. 

b) SED 

The SED points at 850, 450 and 350 pm, shown in Fig. H] 
were obtained integrating the continuum emission over the en- 
velope. We attributed an uncertainty of ~ 30 % to them to ac- 
count for the uncertainty in the envelope size. We also consid- 
ered the IRAS fluxes at 60 and 100 yum, respectively, extracted 
from the IRAS maps at these wavelengths. The evaluation of the 
fluxes was done using the method previously employed for the 
maps at 850, 450 and 350 yum, namely excluding the same re- 
gions (dashed regions in Fig. [TJ to limit the contamination by 
FIR3 and FIR5 and integrating over the rings. We also sub- 
tracted the cloud contribution, estimated to be ~ 0.06 and ~ 
0.07 Jy arcsec"^ at 60 and 100 pm respectively. To account for 
the possible contamination of FIR3 and FIR5 due to the large 
beam of IRAS and the non-sphericity of the source, we took an 
uncertainty of 50 % on the fluxes. Finally, we also considered 
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Fig. 1. Continuum emission maps around OMC2-FIR4 at 850 yum (left panel), 450 jim (middle panel) and 350 jim (right panel). The 
contours mark the continuum flux from 5 % to 75 % of the peak emission by steps of 5 %. The hatched regions have been excluded 
when computing the brightness profile of the FIR4 envelope (see text). The position of the three protostars in the regions, FIR3, 
FIR4 and FIR 5 are marked in the central panel figure. 



the integrated flux at 24 /vm extracted from the Spitzer Space 
Telescope's Multiband Imaging Photometer (MIPS) maps. To 
this end, we retrieved the observations from the Spitzer Science 
archive {http://ssc.spitzer.caltech.edu/archanaly/). The observa- 
tions were obtained the 6th October 2006 as part of the Program 
"Infrared Properties of Edge-on Young Stellar Object Disks" 
(AOR: 30765, PI: Karl Stapelfeldt). The data reduction was per- 
formed using the pipeline S16.0.1. The flux, (5.0 + 2.5 Jy), in 
Fig.|2]was obtained by integration over a 15" radius. 



2.2. Continuum emission: modeling 

To derive the dust physical structure, namely the dust temper- 
ature and dens ity profiles, we used the ID radiative transfer 
code DUSTY dlvezic & Elitzudll997h . Briefly, giving as input 
the temperature of the central object and a dust density profile, 
DUSTY computes self-consistently the dust temperature profile 
and the dust emission. The comparison between the computed 
350, 450, 850 pm brightness profiles (namely the brightness ver- 
sus the distance from the center of the envelope) and SED with 
the observed profiles and SED (see previous paragraph) makes 
it possible to constrain the density profile and, consequently, the 
temperature profile of the envelope. 

To be compared with the observations, the theoretical emis- 
sion is convolved with the beam pattern of the telescope. 
Following the recommendations for the relevant telescope, the 
beam is assumed to be a combination of gaussian curves: at 850 
pm, we use HPBWs of 14.5", 60", and 120", with amplitudes 
of 0.976, 0.022, and 0.002 respectively; at 450 pm, the HPBWs 
are 8", 30", and 12( )" with ampUfiide ratios of 0.934, 0.06, and 
0.006, respectively dSandeU & Weintraubll200lb : at 350 pm, we 



use HPBWs of 12" and 22", w ith amplitude ratios of 0.7, 0.3, 
respectively dHunter et al.ll 19961) . 

We assumed that the envelope density foUows a power law: 

n(r) = no X (1) 

where the power law index, a, is a free parameter of the model, 
as well as the density «o, the density at tq. Besides, the envelope 
starts at a radius R,„ and extends up to R„„,. Both R,„ and Ro„, 
are free parameters of the model. The last input to DUSTY is the 
temperature of the central source, T*, here assumed to be 5000 
K. We verified that the choice of this last parameter does not in- 
fluence the results. Finally, the opacity of the dust as function 
of the wavelength is another parame ter of DUSTY. Following 
nume rous previous studie s ([van der Ta k et al.|ll999HEvans et al.l 
200 Ij; IShirlev et al .' 2002; Young et al. 2003), we adopted the 
dust opacity calculated by Ossenkopf & Henning (1994), specif- 
ically their OH5 dust model, which refers to grains coated by ice. 

In summary, the output of DUSTY depends on a, no, R,„ 
and Ro,„. In practice, the DUSTY input parameters are the 
power law index, a, the optical thickness at 100 pm, tioo, the 
ratio between the inner and outer radius, Y (=Roi/f/Rm) and 
the temperature at the inner radius T,„. The optical thickness 
is, in turn, proportional to the dust column density which 
depends on «() and the physical thickness of the envelope. Note 
that, since the beam of the available maps are relatively large 
(> 7.5" which coiTesponds to a linear length of > 3300 AU), 
the inner region of the envelope is relatively unconstrained 
by the available observational data. In practice, we obtain 
a lower limit to T,„ of 300 K: any l arger value would giv e 
similar results. Finally, as explained in llvezic & Ehtzuri ( Il997h . 
DUSTY gives scaleless results (which makes it very powerful 



4 



N.Crimier et al.: the physical structure of OMC2-FIR4 



because the same grid of models can be applied to different 
sources). This means that to compare the DUSTY output with 
actual observations, it is necessary to scale the output by the 
source bolometric luminosity Lb,,! and the distance. Note that 
the bolometric luminosity is in fact estimated by integrating 
the emission over the full spectrum. By definition, this can 
only be done when the entire SED is known. This is exactly 
one of the outputs of the modeling. So we re-evaluated the 
luminosity of FIR4 iteratively from the best fit model, by mini- 
mizing the We anticipate here that the new value is 1000 
rather than 400 Lq, where we assumed th e most recent esti - 
mation of the distance, namely (437 + 19) pc (iHirota et alj2007l) . 

We run a grid of models to cover the parameter space as re- 
ported in Table [T] The same grid of models were run for four 
values of the illuminating FUV field : Go = 1, 10, 100 and 1000. 
In all cases, w e used the Inter-St ellar Radiation Field (ISRF) 
constructed by lEvans et aTl (|2001[) : combinatio n of the radia - 
tion field introduced bv lBlackl (Il994l) with that of lDrainel(ll978l) . 
Note that, since DUSTY makes the assumption of isotropic scat- 
tering, the computed MIR emission is largely overestimated 
in presence of strong external fields (Elitzur, private commu- 
nic ation). To solve this p roblem, we followed the suggestion 
by lYoung & EvansI (|2005|) to neglect the scattering, artificially 
putting it to zero. The best fit model has been found minimizing 



Parameter 


Range 


a 


0.2-3.9 


Y 


100-2200 


TlOO 


0.1-4.6 


T,„ 


300 K 


T. 


5000 K 



Table 1. Range of the input parameters to DUSTY covered in 
the present study. The range of the a, Y and tioo parameters is 
covered by increasing by 20% their respective value at each step 
of the grid. Note that T,„ and T, are kept fixed as they do not 
influence the results (see text). 



the;i^^ with an iterated two-steps procedure. First, we use the ob- 
served brightness profiles at 350, 450 and 850 i^m to constrain 
a and Y, assuming a value for tioq. Second, we constrain the 
optical thickness tiqo by comparing the computed and observed 
SED, assuming the a and Y of the previous step. The new tioo 
is used for a new iteration and so on. In practice, the iteration 
converges in two steps. This is because the normalized bright- 
ness profiles very weakly depend on tioo, while they very much 
depend on the sizes o f the envelope and on th e slo pe of the den- 
sity pr ofile (see also lJ0rgensen et alj (l2002h and ISchoier et alj 
( I2OO 2')). On the contrary, the optical thickness depends mostly 
on the absolute column density of the envelope, constrained by 
the SED. 



2.3. Results 

We run four grids of models, as discussed separately below: 
a) with a standard illumination FUV field (Go=l) and b) with 
a 10,100,1000 times enhanced field (Go=10,100,1000) (see 
Introduction). In paragraph c), we also discuss why larger G„ 
were not considered, and in paragraph d) we summarize the 
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Fig.4.xlg^i versus Tioo- In these computations, Y is equal to 120 
and a is equal to 0.6. 



results. 
a) Go-1 

Table |2]presents the set of parameters a, Y and tioo, which bet- 
ter reproduce the observations assuming Go=l. Figure |2] shows 
the relevant derived brightness profiles and SED against the ob- 
served ones. Figure [3] shows \ht contours plots obtained by 
considering separately the brightness profiles at 350, 450 and 
850 jUOT, and by combining the three profiles together Figure |4] 
shows the x^ dependence on the rioo parameter 



Observation 


a 


Y 


Tioo 


^red 


V 


850 yum profile 


1.4 


160 




0.72 


10 


450 yum profile 


0.6 


120 




0.63 


10 


350 ytim profile 


0.5 


170 




0.47 


10 


All profiles 


0.6 


120 




1.24 


36 


SED 






0.6 


0.55 


3 



' Note that, in the case of OMC2-FIR4, integrating the model SED 
gives the same results than integrating under the observed SED. 



Table 2. Best fit parameters for the case Go=l. Note Ihaix^^^^ - 
X^/v where v is the number of degrees of freedom. The first line 
reports the best fit obtained using only the 850 fim brightness 
profile; second line, using the 450 jum brightness profile; third 
line, using the 350 fim brightness profile; fourth line gives the 
best fit using the three profiles; the last line gives the best fit 
using the SED. 



The three and;^^^^^ contour plots point to a value of 

Y around 100-200. Conversely, the;^^^^^ and;i'45g contour plots 
constraint a to a lower value than 1, around 0.5-0.6, whereas 
the Xgjo would rather indicate a larger value for a, although 
the value 0.6 is sti l l ac ceptable. Note that the solution found 
by 'J 0rgensen et al.l (1200 6) relies on the 850 fim profile only, 
and, therefore, gives a large a value, consistent with our xI^q 
plot. The XsED Pl^*^ (FigHJ points to a value of tioo of 0.6. In 
minimizing the x\ed' varied the source luminosity from 
400 to 1500 Lq. The best fit is obtained for a source luminosity 
equal to 1000 Lq. 
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Fig. 2. Observed brightness profiles at 350 (upper left panel), 450 (upper right panel), 850 ;um (lower left panel) and SED (lower 
right panel). The curves report the best fit obtained in the two cases Go=l (solid line) and 1000 (dashed line). The dashed-dotted 
lines represent the beam pattern of the telescope adopted at 350 450 and 850 fim. Note that the SED plot reports the ISO-LWS 
spectrum between 45 and 200 fim for completeness, although it has not been considered in the;^^^ analysis, due to the relative larger 
caUbration uncertainty compared to the IRAS data. 



b) Go^lO,100,1000 

The best fit values of or and Y for cases of an enhanced illu- 
mination UV field are presented in Figure |5] The first thing to 
notice is that the xlu does not change appreciably for G,, equal 
to 1,10,100 or 1000: the minimum value is 1.24, 1.23, 

1.19 and 1.20 for Go=l, 10, 100 and 1000 respectively. In other 
words, the available continuum observations, both the profiles 
and the SED, cannot distinguish which of the four models is 
better. Furthermore, Figure |5] shows that the x'au contour plots 
point to the same Y and a values. Similarly, the tioo value is 
0.6 for the four cases Go=l, 10, 100 and 1000. The situation is 
illustrated in Fig. |2] where the best fit predictions are compared 
to the observations for the two cases G„=l and 1000. Both 
models reproduce fairly well the observations, as implicit in the 
best-;^'^ similar values. Note, however, that the Go=1000 case 
predicts slightly larger fluxes, due to the enhanced temperature 
at the border of the envelope. 

c) larger Go 

We did not explore in detail the case of larger Go for three 
reasons. The first one is that previous line observations showed 
that the FUV field in the O MC2 region is " o nly" 5 00 times 
the Interstellar field. Indeed, [Herrmann et al.l J 19971) mapped 
the OMC-2 cloud in the CII-157 OI-63 and -146 yum 
lines with the spectrometer FIFI on board the Kuiper Airborne 



Observatory. They detected extended emission associated with 
the Photo-Dissociation Region (PDR) enveloping the whole 
OMC-2 molecular cloud. These authors concluded that OMC-2 
is illuminated by a FUV field whose intensity is Go ~ 500. Note 
that this is the FUV field impinging on the cloud and that the 
eff'ective Go seen by the FIR4 envelope is probably lower than 
this. The second reason is that varying Go from 1 to 10"^ does not 
improve the;^'^ value. The third reason is that the Go =10^ case 
suffers of severe convergence problems, and it was not possible 
to derive a large enough number of runs for a meaningful 
analysis. 

d) Summary of the adopted solution 

Table|3]summarizes the value of the best fit parameters, obtained 
by considering all the profiles and the SED contour plots and 
assuming the Go=l case. Some relevant physical quantities are 
quoted in the same table. Fig.|6]shows the dust density and tem- 
perature profiles of the best fit models with Go=l and 1000 re- 
spectively. Note that the dust temperature in the skin of the enve- 
lope is larger by ~ 20-30 K in the case Go=1000 with respect to 
the Go=l case. This incr ease concerns a r elative ly small region, 
of a few thousand AU. |j0rgensen et alj (l2006l) found a larger 
warm region, of about lO'* AU, because of the steeper adopted 
density distribution {a-2): in this case, the FUV photons can 
penetrate deeper into the envelope. 
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Fig.3.xjgd contour plots iY,a) for the models with Gg-l- In these computations, tioo is equal to 0.6. The contours show the loci 
of Xhe. x^red Values equal to 1.1, 1.5, 2.5 and 5 times the minimum The upper left panel is obtained by comparing the model 
predictions with the 850 yum brightness profile; the upper right panel refers to the 450 jum profile; the lower left panel refers to the 
350 jum profile; the lower right panel makes use of the three profiles. 



Fixed input parameters 



Distance, d 


437 pc 


Stellar temperature, T^, 


5000 K 


Dust temperature at r,, Ti{r\) 


300 K 


Dust opacity (0H5) at 100 yum, kiqq 


86.5 cm^g-' 


Best fit parameters 


Luminosity, L 


1000 Lq 


Dust optical depth at 100 /im, Tjoo 


0.6 


Density power law index, a 


0.6 


Envelope thickness, rom/r; 


120 


Physical quantities 


Inner envelope radius, Hn 


lOOAU 


Outer envelope radius, 


12000AU 


Radius at T^,,,, = 100 K, riooK 


440 AU 


H2 density at riooK, «o 


4.3 X 10« cm-3 


Envelope mass, Menv 


30 M© 



Table 3. Summary of the dust radiative transfer analysis of 
OMC2-FIR4. The first part lists the fixed input parameters, the 
second part reports the best fit parameters, and some relevant 
physical quantities corresponding to the best fit model are re- 
ported in the third part. 



3. Gas temperature profile 

3. 1 . Model description 

rCeccarel li et al.l(ll996t) . lDoty & Neufeldl(ll997h and lMaret et aP 
(12002) showed that the gas is thermally decoupled from dust in 
the inner regions of low and high mass protostellar envelopes. 
The reason for that is the large water abundance in the gas phase 
caused by the sublimation of the grain mantles. The same phe- 
nomenon may occur in the envelopes of intermediate mass pro- 
tostars, so we explicitly computed the gas temperature profile of 
the envelope surrounding FIR4. For that we explicitly computed 
the equilibrium temperature by equating the gas cooling and 
he ating terms at e a ch rad ius. Following the method described 
in ICeccarelli et al.l (Il996l) . we considered heating from the gas 
compression (due to the collapse), dust-gas collisions and photo- 
pumping of H2O and CO molecules by the IR photons emitted 
by the warm dust close to the centefl The cooling is mainly 
due to rotational lines from H2O and CO, plus the fine struc- 
ture lines from O. Therefore, the gas temperature depends on the 
abundance of these three species. In practice, though, only the 
water abundance is a real parameter of the model, because the 
CO and O lines are optically thick in the range of the CO and 
O abundances typical of protostellar envelopes. For this reason. 



^ Cosmic rays ionization is a minor heating term in the protostellar 
envelopes. 
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Fig.5. contour plots {Y,a) for the models with Go=l (black lines), 10 (blue lines), 100 (green lines) and 1000 (red lines). In 
these computations, tioo is equal to 0.6. The contours show the loci of the x^^j values equal to 2.5 times the minimum From 
top to bottom: contours of the 850 fim, 450 yum, 350 yum and the three together. 



we computed various cases for the water abundance, as it is gen- 
erally poorly constrained in protostellar envelopes, and totally 
unconstrained in FIR4 (see ® . We adopted a step function for 
the water abundance profile to simulate the jump caused by the 
ices sublimation. The jump is assumed to occur at 100 K. We 
considered the H2O abundance (with respect to H2) X(H20)o,« 
in the outer envelope, where T< 100 K, equal to 10"^, 10"^ and 
10"^. We also considered three cases for the abundance in the 
inner region X(H20)„,, 10"^, 10"^ and 10"^. Finally, we studied 
the case with Go = 1000. The run parameters are summarized in 
Table H 



Model 


X(CO) 


X(0) 


X(H20)„„, 


X(H20),„ 


Go 


r 


1 X 10-* 


5 X 10"* 


1 X 10"' 


1 X 10"' 




2 


1 X lO"'* 


5 X 10"* 


1 X 10"* 


1 X 10"' 




3 


1 X 10"'* 


5 X 10"* 


1 X 10"* 


1 X 10"* 




4 


1 X 10"** 


5 X 10"* 


1 X 10"' 


1 X 10"' 




5 


1 X 10"** 


5 X 10"* 


1 X 10"' 


1 X 10"' 




6 


1 X lO"'* 


5 X 10"* 


1 X 10"' 


1 X 10"* 




7 


1 X 10"** 


5 X 10"* 


1 X 10"* 


1 X 10"' 


1000 



Table 4. The different run models (column 1). Column 2 to 5 
report the adopted abundances of the main gas coolants: X(CO), 
X(0) and the H2O abundance in the T>100 K region X(H20),„ 
and outer region X(H20)o«/. Column 5 reports the FUV illumi- 
nating field Gq. Note: " Model 1 is the reference for the studies 
of the water line spectrum presented in ^ 



To co mpute the cooling froni the lin es we used the code de - 
scribed in lCeccarelli etalld 19961 llool and lParise etan ( l2005h . 
The same code has been used in several past studies, whose re- 
sults have been substantially confirmed by other groiips (e.g. the 
analysis on IRAS 16293-2422 by (Schoier et al.ll200a) ). Briefly, 
the code is based on th e escape probabi l ity for malism in pres- 
ence of warm dust (see iTakahashi et al.l ( Il983l) ). where the es- 
cape probability /3 is computed at each point by integrating the 
line and dust absorption over the solid angle Q as follows: 



kh + kd (^L + ^d) 



r 1 -exp[-(^L + ^d)ALth] 
2 I ^J" ~ (2) 



ALd 



where k]^ and kd are the line and dust absorption coefficients re- 
spectively, and ALth is the line trapping region, given by the fol- 
lowing expressions: 



ALd 



2Av,h - 



1-V 
2^ 



(3) 



in the infalling region of the envelope (where arcos(//) is the 
angle with the radial outward direction) and 



ALd 



(4) 



in the static region (where /?env is the envelope radius). In the 
present calculations, we assumed that the entire envelope is col- 
lapsing in free-fall towards a central object of 2 Mq. In practice. 



8 



N.Crimier et al.: the physical structure of OMC2-FIR4 




10^ 7 

I 

e 




100 1000 10000 

Radius (AU) 

Fig. 6. Dust temperature (upper panel) and H2 density (lower 
panel) profiles from the best fit obtained in the two cases Go=l 
and 1000. The plain line and the dotted line represent the cases 
Go=l and 1000 respectively. 



the photons emitted by the dust can be absorbed by the gas and 
can pump the levels of the water molecules. This, indeed, is an 
important factor in the population of the water levels, and, for 
the highest energy levels, even the dominant one (®. In addi- 
tion, H2O and CO molecules can be pumped by absorption of 
the NIR photons emitted by the innermost warm dust. Since the 
densities and temperatures of the regions of the envelope tar- 
geted by this study are not enough to populate the levels at the 
vibrational states, the effect of the NIR photons is an extra heat- 
ing of the gas, as described in the Ceccarelli et al. (1996) article. 
Note that the code takes into account the dust with temperatures 
up to 1500 K, by following the algorithm described in Ceccarelli 
etal. (1996). 

For the collisional coeffic ients of w a ter w ith hydrogen 
molecules, we used the data bv lFaure et al.l ( l2007h available for 
the temperature range 20-2000K. This data set includes quasi- 
classical results for the highest rates (those larger than 10"'^ 
cm^s ' ) and quantum scaled H20-He results for the lowest rates. 
Recent quantum calculations on ortho-H20 by Dubernet and co- 
workers have shown that the quasi-classical rates can be in error 
by as much as a factor of 100 but that, in general, they are accu- 
rate to within a factor of 1-3 (Dub ernet et al. 2009). It should be 
noted that the rates of iFaure et al.l (f2007l) are currently the only 
complete and consistent set of data for both ortho- and para- 



Fig. 7. The gas temperature profile of the collapsing envelope 
of OMC2-FIR4. The different curves refer to different values of 
the inner envelope water abundance X(H20)„,: 1 x 10"^ (dotted), 
1 X 10"^ (dashed) and 1 x 10"^ (dotted-dashed) respectively. In 
these computations, X(H20)om/ is 1 x 10"^ and Go=l. The solid 
line refers to the dust temperature profile. 

H2O colliding with H2. We also note that these rates have been 
recently extrapolate d in order to cover energ y levels and temper- 
atures up to 5000K jFaure & Josselin|[2008 ). Since the ortho to 
para conversion process of H2 is chemical rather than radiative, 
the ortho-to-para ratio H2 OPR is highly uncertain in the inter- 
stellar medium. Here we assume that in warmer gas it is in Local 
Thermal Equilibrium and, therefore, follows the Boltzmann dis- 
tribution: 

(24 + 1)1(2/ +l)exp(-^) 

OfK - ^-7- (5) 

(2/^ + 1)1(27+ l)exp(-^) 

where and Ip are the total nuclear spin, corresponding to 
whether the hydrogen nuclear spins are parallel (/„ = 1, ft) or 
anti-parallel (Ip = 0, ti)- The sum in the numerator and de- 
nominator extends over all ortho and para levels J, respectively. 
Similarly to H2, water comes in the ortho and para forms. In this 
case, since the water is the dominant gas coolant only in the re- 
gions where the dust temperature exceeds 100 K, we assumed 
OPR equal to 3, strictly valid for gas temperatures larger than 60 
K. Since the water lines are optically thick, the cooling depends 
on the velocity field, assumed to be that of an envelope collaps- 
ing in free-fall towards a central object of 2 M0 (see above). 
We checked the influence of our results against this assumption, 
running a case with a constant velocity field of 0.5 km/s. The 
difference in the gas temperature between the two cases never 
exceeds 10 K. 

3.2. Results 

Figure [7] shows the computed gas temperature profile obtained 
with different values of X(H20)„, in the case Go = 1. Figure 
[8] shows the different contributions to the heating and cooling 
rates. Similarly t o what had been found in low mass protostars 
jCeccarelli et alj[r9 96, 2000: lMaret et al.ll2002l) . the gas temper- 
ature tracks the dust temperature in the outer envelope while gas 
and dust are decoupled in the inner part of envelope, where the 
icy grain mantles sublimate. The heating is dominated by com- 
pression of the collapsing gas across the entire envelope, even 
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Fig. 8. Heating (top panel) and cooling (bottom panel) rates as 
function of the radius, computed assuming that the inner H2O 
abundance is equal to 1 x 10"^ while the outer abundance is 1 x 

io-«. 



though the dust-gas collisions becomes comparable to the com- 
pression heating in the inner envelope. Although important in 
the very inner regions, the H2O photo-pumping never dominates 
the heating contrarely to what happens in the studied low mass 
protostars. The cooling, on the other hand, is dominated by H2O 
line emission in the inner envelope, by the 01 line emission in 
the intermediate region and by the CO in the outermost regions 
of the envelope. Note that the increased water abundance causes 
an increased cooling of the gas, which brings the equilibrium 
gas temperature to lower values than the dust temperature. This 
phenomenon, already predicted in low mass protostars, is much 
more marked in the FIR4 case, leading to more than 100% of 
difference (with respect to the gas temperature) in the dust and 
gas temperatures for the case of the highest water abundance 
(1 X lO"'*). For example, at 100 AU the dust temperature is 300 
K, whereas the gas temperature varies from 200 to 80 K depend- 
ing on the assumed X(H20),„, 1 x 10"^ and 1 x lO"'* respectively. 
The phenomenon is more marked in FIR4 than in the studied low 
mass protostars because of the relatively lower density of the re- 
gion where the icy grain mantles sublimate in FIR 4 than in the 
low mass protostars, or, in other words, because the FIR4 enve- 
lope is warmer. Note that we obtain similar results also for larger 
illuminating FUV fields. 



Fig. 9. Synthetic rotational diagram derived from the water line 
emission using the reference model (Model 1, Table |4| inte- 
grated over the whole envelope. Crosses and diamonds trace the 
ortho and para water, respectively. 

We emphasize that this result is a consequence of the derived 
shallow dependence of the density distribution, which is con- 
strained from the fit of the maps. The dependence is strictly valid 
only at scales larger than the smaller telescope beam, namely 8" 
(equivalent to a radius of about 1700 AU) and the SED fit only 
gives the total column density, which, coupled with the density 
dependence on the radius (constrained by the maps), constrains 
the density at these scales. While we cannot exclude the presence 
of a denser compact object hidden by the envelope, it seems un- 
likely that the envelope density gradient increases inwards, be- 
cause this would be unphysical. 

Clearly, the water abundance in the inner region of FIR4 will 
have a great impact not only on the emerging water spectrum 
but also on the emerging line spectrum of any molecule (abun- 
dant in the inner region), and has to be correctly taken into ac- 
count to give reliable molecular abundances. Conversely, given 
the large effect, in principle appropriate multiline observations 
of any molecule will be able to constrain the inner region water 
abundance and the present model predictions. Note that varying 
the outer abundance X(H20)o«f does not have effect on the gas 
temperature, as in the outer region the cooling is dominated by 
the CO and O lines. 

4. Predicted water line spectrum 

4. 1 . Reference model 

Here we report and discuss the predicted spectrum of our ref- 
erence model. Next paragraph will discuss how it depends on 
the parameters of the model. We adopted the Model 1 of Table 
ID as reference model . We first discuss the general water line 
spectrum by means of the synthetic rotational diagram, and then 
we discuss the specific predictions for the two spectrometers on 
board Herschel: HIFI and PACS. 

Figure|9]shows the synthetic rotational diagram derived from 
the line emission integrated over the whole envelope. As ex- 
pected, the theoretical points do not lie on a compact and straight 
line, reflecting the different line optical depths, the gradients in 
density and temperature of the envelope and non-LTE effects. An 
illuminating example is represented by the fundamental transi- 
tions of the ortho and para water lines at 557 and 1113 GHz 
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respectively. We will discuss these two lines in detail because, 
first, they will certainly be important observational diagnostics 
and, second, they offer a great pedagogic case. The situation is 
illustrated in Fig[TOl where we report the profile of the emission 
of the ortho and para Hj*"© and Hj^^O fundamental lines as func- 
tion of the radius. Figure[TT] with the beta escape probability as 
function of the radius for the two fundamental Hj^O lines, also 
greatly helps to interpret the emerging line fluxes for the two 
lines. The ortho-Hj^O fundamental line emission (FigfTOb peaks 
at the border of the envelope and it decreases inwards because 
of the decreasing emitting volume. The para-H20 fundamental 
line shows approximately the same behavior. If the lines were 
optically thin and LTE populated, the expected flux ratio of the 
para over ortho fundamental line would be between 3 and 4 for 
a temperature between 50 and 200 K. Any departure from this 
value originates from a combination of line opacity and non-LTE 
effects. In the outer region the ratio is lower than 1 : the para-H20 
line is optically thin, whereas the ortho-H20 lines is moderately 
optically thick (FigfTTTi. Therefore, the much lower emission of 
the para-HiO line with respect to the ortho-line is due to the 
non-LTE population effect, more accentuated in the para-H20 
line. The situation is reversed in the inner region, where ices sub- 
limate: the para-H20 fundamental line becomes about ten times 
brighter than the ortho-H20 fundamental line because of the line 
opacity, which is much larger in the ortho-H20 line than in the 
para-H20 line (FigfTTTi. In fact, the increase in the water abun- 
dance by a factor 1000 gives rise to a jump in the line emission 
by a factor 3 in the ortho-H20 line and 30 in the para-H20 line, 
and this can only be due to the larger opacity of the ortho-H20 
line as the excitation conditions do not change when ices sub- 
limate. In summary, the emission from the water lines is due, 
in principle, to a rather complex combination of line opacity, 
non-LTE effects and emitting volume (namely temperature and 
density gradient). Evidently, the intensity ratio of lines from the 
Hj^O and Hj^O isotopologues is far to give the "opacity" of the 
line, as it is a combination of the penetration of the line and the 
opacity itself. 

Table|5]lists the predicted water line fluxes for the two spec- 
trometers on board Herschel: HlFl and PACS. Note that, in both 
cases, we computed the signal after convolving the theoretical 
line intensity map with the instrument beam which vary from 
39" to 13" with the frequency varying from 500 GHz to 2000 
GHz (HIFl frequency range) and from 13" to 5" for wavelengths 
from 210 jum and 60 fj.m (PACS wavelength range). 



Based on the (preliminary) sensitivities reported 
on the Herschel Observation Planning Tool HSpot: 
http://herschel.esac.esa.int/Tools.shtml), several ortho and 
para lines are predicted to be detectable by the two Herschel 
spectrometers: about a dozen in the HIFI frequency range and 
twice more in the PACS wavelength range. The Hj**© ortho and 
para lines are also predicted to be detectable by HlFl, and 100 
and 20 times less bright than the respective lines of the Hj^O 
respectively. Note the counter-intuitive result: the para-Hj^O 
line seems to be more optically thick than the ortho-Hj^O line! 
As explained above, this is not the case, of course: the line 
intensity ratio (from which the line optical depth is usually de- 
rived) is due to the combination of optical depth plus excitation 
(non-LTE) effects, and the final result is not easily predictable. 
In our reference model, no observable line is predicted to be in 
absorption. 




1000 10000 
Radius (AU) 




t3 
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Fig. 10. Emission profile, of the ortho (solid line) and para 
(dashed line) water fundamental lines at 557 and 1113 GHz re- 
spectively as function of the radius, in the case of the reference 
model (see Table HJi. The plotted quantity is the contribu- 
tion of each shell to the flux integrated over the whole envelope. 
H2^0 and H^^O emission profiles are represented on the top and 
bottom panel, respectively. 



4.2. Other models 

Here we explore the sensitivity of the results reported in the pre- 
vious paragraph against the variation of the three main parame- 
ters of the model: the water abundance in the inner (X(H20),„) 
and outer (X(H20)o„,) envelope, and the illuminating FUV field 
Go. 

Figures [12] and [13] show the ratio between the line intensities 
of the reference model (Model 1 of TableHJ and the line intensi- 
ties predicted by models with different X (H20)out and X(H20)if i 
respectively. As n oted by other authors (ICeccarelli et al.l [20001 : 
iMaret et ail l2002t) . lines with upper level energies lower than 
about 200 cm"' are sensitive to X(H20)o«/ and insensitive to 
X(H20),„, because these lines mostly originate in the outer en- 
velope for excitation and line opacity reasons. A variation of a 
factor 10 in X(H20)o,„ leads to an almost similar variation in 
the line intensity of the lowest lying lines. The higher the up- 
per level energy the smaller the variation. Conversely, lines with 
upper level energies larger than about 200 cm"' are sensitive to 
X(H20),„ and insensitive to X(H20)o«r. In this case, variations 
by a factor 10 in X(H20),„, going from 1 x lO"'' to 1 x 10"^ 
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Fig. 11. Ratio of the H\^0 escape probability over the H^^O 
escape probability of the ortho (solid line) and para (dashed 
line) water fundamental lines obtained with the reference model 
(model 1 of Table |4| at 557 and 1113 GHz respectively as func- 
tion of the radius. 



can lead to variations in the lines fluxes even 100 times larger. 
This extreme variation, 10 times larger than the difference in the 
X(H20),>, variation, occurs to some lines in the 50-200 fim wave- 
length range. This phenomenon is due to the fact that those lines 
are in absorption rather than in emission in the region just after 
the ices sublimation, resulting in an additional decrease of the 
emerging line flux. The higher the X(H20),„ the smaller the ab- 
sorption depth. When X(H20),„ reaches 1 x 10"^ the absorption 
region generally vanishes. In addition, many high lying lines are 
prevalently populated by absorption of the photons emitted by 
the dust, so that they are particularly sensitive to the dust contin- 
uum. 

Note that as mentioned in the Sect. 12.21 the inner region of 
the envelope is relatively unconstrained by the available obser- 
vational data. Therefore we derived water line spectrum predic- 
tions varying the power law index of the density profile a of 
about ~ 30% in the inner part. We observed a variation of the 
line intensity of a factor 5-10 for the transitions with upper 
level energy > 300-400 cm ' and lesser than 2 for the lower 
lines. Finally, the predicted line intensities do not vary apprecia- 
bly when the illuminating FUV field changes from 1 to 1000. 
Therefore, observations of water lines will be extremely helpful 
in constraining the water abundance across the envelope, but will 
not be sensitive to the illuminating FUV field. 

4.3. Effect of gas-dust thermal decoupling 

As presented in ^ the large quantity of water vapor injected 
into the gas in the inner part of the envelope causes a dramatic 
decoupling between the dust and gas temperatures (see Fig. |7]i- 
Obviously, this effect has a great impact in the interpretation of 
the water line emission. This is illustrated in Figures [14] where 
we report the ratio of the water line intensities obtained by con- 
sidering the gas temperature self-consistently computed (model 
1) over the case where Tj,,, is assumed to be equal to Tdusi- 
Assuming artificially Tgas = Tdust leads to differences in line 
fluxes up to two orders of magnitude. Since the decoupling oc- 
curs in the inner part of the envelope, the larger the upper level 
energy the larger the difference, for both ortho and para lines. 
Note that the fluxes of the two fundamental ortho and para lines 



PACS range 


Transition 


Wavelength 


Flux 




(pm) 


(erg s"' cm" 2) 








221 — > 330 


66.44 


2.99E-14 


220 — * 331 


67.09 


1.13E-14 


3o3 —'330 


67.27 


2.28E-14 


2l2 321 


75.38 


1.55E-13 


3l2 423 


78.74 


4.51E-14 


2ll — > 322 


89.99 


5.09E-14 


4o4 5|5 


95.63 


1.05E-14 


4i4 5o5 


99.49 


4.21E-14 


111 — > 220 


100.98 


7.80E-14 


llO — * 221 


108.07 


1.43E-13 


3o3 — * 4i4 


113.54 


8.45E-14 


423 -> 432 


121.72 


1.14E-14 


3l3 -> 4o4 


125.36 


5.34E-14 


4l4 -» 423 


1 QO /I 1 


/.yzrl-i4 


321 — » 330 


136.49 


1.35E-14 


2o2 — » 3l3 


138.53 


7.84E-14 


3l3 — > 322 


156.20 


2.25E-14 


2l2 — > 3o3 


174.62 


7.42E-14 


loi — > 2i2 


179.53 


1.17E-13 


2l2 — > 221 


180.49 


4.97E-14 


HIFI range 


Transition 


Frequency 


Flux 




(GHz) 


(KKms-') 


HfO 






loi — > lio 


556.96 


14.8 


2o2 — > 2ii 


752.04 


1.06E+00 


In — » 2o2 


987.95 


1.31E+00 


3o3 — » 3l2 


1097.34 


1.46E+00 


Ooo — > 111 


1113.35 


4.32E+00 


221 — * 3i2 


1153.09 


2.50E+00 


3l2 — * 321 


1162.93 


6.97E-01 


2ii — » 220 


1228.81 


5.33E-01 


4l3 -» 422 


1207.62 


2.45E-01 


5l4 ^ 523 


1410.65 


1.47E-01 


2l2 221 


1660.99 


2.33E+00 


4o4 -» 4,3 


1602.23 


1.46E-01 


loi — > 2|2 


1669.87 


5.46E+00 


2l2 — > 3o3 


1716.83 


3.37E+00 


523 — ' 532 


1867.75 


3.13E-02 








loi — » lio 


556.96 


1.88E-01 


Ooo — » 111 


1113.35 


1.48E-01 



Table 5. Predictions of the line fluxes (after subtraction of the 
continuum) of the water lines observable with the Herschel spec- 
trometers, HIFI and PACS. The predictions refer to the reference 
model (model 1 of Table |4]i. 



are not affected by the Tgas - Tdust choice. Finally, we did the 
same study with the model 2 ( X(H20),„ = lO"^). In this case, 
due to the smaller X(H20),„, the decoupling is less important 
than in the previous one, leading to changes in lines intensities 
up to one order of magnitude. 

4.4. Constraints from ISO data 

Observations of FIR4 were obtained by the spectrometer ISO- 
LWS in the grating mode (spectral resolution ~200) and Fabry- 
Perot (spectral resolution ~ 10"*). We retrieved the data from 
the ISO Data Archive (http://iso.esac.esa.int/ida/). No water lines 
are detected. Table |6] summarizes the upper limits obtained for 
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ivioQei 1 


ivioaei 


TQM 




Wavelength 


Transition 


Intensity 


Intensity 


Instrument 


Upper Limit (3cr) 


(pim) 




(erg.s '.cm ^) 


(erg.s ^ .cm ^) 




(erg.s '.cm ~) 


ortho 












75.4957 


854 — > 845 


2.6E-20 


1.9E-19 


LWS04 


9.78e-12 


108.073 


221 ^ lio 


1.3E-12 


1.5E-12 


T \T TC^ (\ A 

LWS04 


2.21e-12 


1 7Q ^^97 
1 ly.JZ. 1 


'-11 ^ ioi 




7 OP 1 9 


T WQOzL 


9 ^f\(^ 1 1 
z. joe- 1 1 


113.538 


4l4 3o3 


8.3E-13 


8.7E-13 


LWSOl 


1.2e-ll 


para 












100.983 


220 111 


7.3E-13 


9.1E-13 


LWSOl 


8.37e-12 


138.527 


3l3 — > 2o2 


7.3E-13 


7.3E-13 


LWSOl 


4.14e-12 


156.197 


322 3l3 


2.2E-13 


3.4E-13 


LWSOl 


1.60e-12 



Table 6. The brightest lines predicted by models 1 and 6 (the model with the largest water abundance) compared with the upper 
limits derived by the ISO observations.) 
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Fig. 12. Ratio between the line fluxes of the reference model 
(model 1 of Table |4| and the line fluxes predicted by models 
with same X(H20)i„ but different X(H20)o„f. Flux ratios be- 
tween model 4 (model 5) and model 1 are represented by squares 
(triangles). Filled and empty symbols refer to the lines emitted in 
the PACS and HIFI bands, respectively. The upper (lower) panel 
reports ortho (para) water line intensities ratios. 



Fig. 13. Ratio between the line fluxes of the reference model 
(model 1 of Table HJ and the line fluxes predicted by models 
with same X(H20)o,,/ but different X(H20),„. Flux ratios be- 
tween model 2 (model 3) and model 1 are represented by squares 
(triangles). Filled and empty symbols refer to the lines emitted in 
the PACS and HIFI bands, respectively. The upper (lower) panel 
reports ortho (para) water line flux ratios. 



the predicted brightest lines together with the predictions of the 
reference model and model 6. Unfortunately, the ISO sensitivity 
is not enough to put sensible constraints to the water abundance 
across the FIR4 envelope. 



5. Concluding remarks 

We have analyzed in great detail the continuum emission from 
the Intermediate Mass protostar OMC2-FIR4, with the aim of 
deriving the physical structure of its envelope, a mandatory first 
step for further studies to understand the formation process. Our 
analysis led to a new estimate of the FIR4 luminosity, which 
is 1000 Lq. The density of the envelope surrounding FIR4 has 
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Fig. 14. Ratios between the line fluxes of the reference model 
(model 1) and the reference model with gas-dust non-thermally 
decoupled (namely Tgas^Tdust, as function of the upper level en- 
ergy of the transition. Upper and lower panels show ortho and 
para H2O lines, respectively. 



a shallow dependence on the radius, the density power law in- 
dex being only 0.6. Since systematic studies of the IM pro- 
tostars envelopes have not been published yet, we can tenta- 
tively compare the FIR4 envelope with low and high mass pro- 
tostellar env elopes, where similar st udies have been caiTied out. 
Specifically. |j0rgensen et afl ( |2002|) analyzed 18 Class and I 
sources and found that the average power law index a in Class 
sources is 1.3 ± 0.4 while in Class I sources it is 1.7 + 0.1, 
significatively lar ger than the value we found in FIR4. Similarly, 
Ivan der Tak et al.l (|2000) studied a sample of high mass proto- 
stars and found a = 1 .4 + 0.4. One has to notice that, however, 
there are exceptions, with sources in both low and high mass 
showing a smaller than unity a value: LI 527 (a = 0.6) and 
L483 (a = 0.9) in the Class sources, GL7009S (a = 0.5) 
in the high mass protostars sample. It is not clear what makes 
these sources "anomal ous": the presence of strong asymmetries 
(lJ0rpensen et al.ll2002h have been suggested as possible reason. 
The case of FIR4 seems to fall in this "anomalous sources" cat- 
egory, and further studies, possibly on chemistry, are needed to 
say more. 

Giving the suggestion bv lJ0rgensen et al.l (l2006h that a strong 
FUV field (G()=l X lO'*) illuminates the FIR4 envelope, we ex- 
plored the cases of different FUV fields. As already noted by the 
same authors, however, the dust continuum cannot really distin- 
guish whether a strong illuminating FUV field is impinging on 



the envelope. In fact, |j0rgensen et al.l (l2006l) adopted a steeper 
density distribution (a-2) which allows the FUV photons to pen- 
etrate deeper into the envelope. Their conclusions were based 
on submillimeter lines from CO and H2CO, which would be 
exceedingly bright if they were emitted in the envelope. They 
attributed the lines to the warm gas at the border of the enve- 
lope, heated up by the hypothetical larg e FUV field. However , 
as discussed in ^2.31 OI and CII maps bv lHerrmann et al] d 19971) 
showed that the entire OMC2 region is illuminated by a Go=500 
FUV field, which would imply an even lower FUV field on 
the FIR 4 envelope. One has also to notice here that large scale 
maps bv ISchloerb & LorenI (1 1 982' ) show that the CO (1^0) line 
is b right (~40 K) over th e whole OMC2 region, a fact that 
lead [Herrmann et alj ( 1 19971) to attribute the CO emission to the 
PDR associated with the cloud. In addition, several outflows are 
known to "pollute" the CO emission in the region, in partic- 
ular the one originating from FIR3 (25" North of FI R4: Fig. 
[B and reaching FIR4 and FIR5 ( Williams et alj|2003h . All the 
above considerations together lead to conclude that the FUV 
field impinging FIR4 is not anomalously large and less than 500. 
Therefore, giving the presence of a bright PDR and a "polluting" 
outflow from FIR3, caution is needed in interpreting the low ly- 
ing water lines, as much as lines from any molecule, separately 
from the whole molecular cloud emission. 

One major motivation of the present work is the prediction 
of the water line spectrum from FIR4, as this source will 
be observed in the 500-2000 GHz frequency range by the 
incoming Herschel Space Observatory (FIR4 is a target of 
the Key Program "HIFI Spectral Surveys of Star Formation 
Regions": http://www-laog.obs.ujf-grenoble.fr/heberges/hs3f/). 
In the present study, we have shown that water is indeed a 
key molecular species, because of its great impact on the 
gas cooling in the region where the dust temperature exceeds 
100 K, the sublimation temperature of the dust grain ices. 
The large quantity of water vapor injected into the gas by the 
sublimated ices very efficiently cools the gas, causing a dramatic 
decoupling between the dust and gas temperatures. Depending 
on the abundance of the injected water vapor, the difference in 
the temperature can be as high as 50 K at the sublimation radius 
(namely 50%!) and even larger going inward. For example, at 
100 AU the dust temperature is predicted to be around 280 K 
whereas the gas temperature is 80 K if the water abundance is 
1 X lO"'*. Obviously, this has a great impact in the interpretation 
of the water line emission as much as the emission from any 
molecular species emitting in the inner region. In fact, the 
comparison of the water line emission between the case where 
dust and gas are assumed to be thermally coupled and the 
case where the gas temperature is self-consistently computed 
shows that the difference of the line intensity can reach two 
orders of magnitudes for lines with large upper level energies 
(namely the lines excited in the innermost region, where gas 
and dust decouple). Therefore, our important second conclusion 
is that caution has to be applied in interpreting the line emission 
from FIR4, as much as any source with a similar luminosity 
and envelope structure. Gas and dust temperature can be very 
different and in order to derive correct molecular abundances 
(including water abundance) they have to be both estimated, 
accounting for all terms of heating and cooling. Avoiding that 
may lead to very wrong conclusions. 
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